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INTRODUCTION 

The atmospheres of giant and supergiant 
stars of the M-, S-, C-type appear to be intrin- 
sically very different from the late-type dwarfs. 
Although M dwarfs possess a chromosphere, 
a transition region, and a corona that are qual- 
itatively similar to that of the Sun, these lum- 
inous giant and supergiant objects appear to be 
devoid of transition regions and coronae, but 
to be characterized by vast amounts of circum- 
stellar matter (Rieu, Lefevre, and Glassgold and 
Huggins, this volume) and often important 
mass loss (Goldberg, this volume). As shown 
by observations of highly ionized lines in the 
ultraviolet with the International Ultraviolet 
Explorer (IUE), there are no indications that 
plasmas of temperatures of 10 5 K, correspond- 
ing to a transition region, exist. A recent review 
on late-type IUE observations is given by Ba- 
liunas (1983). The same negative results were 
obtained by absence of direct detections of the 
X-ray fluxes by means of the Einstein Obser- 
vatory Satellite, indicating the absence of co- 
ronae (Vaiana et al., 1981). Indirect observa- 
tions of coronae made by observing the He I 
line at 10830 A (Zirin, 1982) also gave negative 
results. This chromospheric line can be formed 
by a photoionization/recombination process, 
produced by illumination of a coronal radia- 


tion. Moreover, no lines from high stages of 
ionization were seen in a short wavelength ex- 
posure with IUE of the carbon star, TX Psc, 
with an exposure of 7 hours (Johnson and 
O’Brien, 1983). 

Some recent evidences show that, above the 
photosphere, a vast chromosphere could exist 
(Stencel, 1982a), at least for the early red giants 
and supergiants, extending to several stellar 
radii. These chromospheres connect the photo- 
sphere to the circumstellar regions. The im- 
portant mass loss observed in these stars begins 
to operate somewhere in these chromospheres. 
If these chromospheres are extended, a new 
problem appears: to find the mechanism that 
could support them— turbulence? magnetic 
field? The answers are not definite. 

In our discussion and construction of quasi- 
thermal models, we will relax the local thermo- 
dynamic equilibrium (LTE) assumption in 
photospheres and chromospheres. In general, 
hydrostatic equilibrium will be valid 
everywhere. Radiative equilibrium is valid in 
the photosphere, but this will no longer be the 
case in the chromosphere. Some velocity fields 
are, in this context, sometimes “ad-hoc” 
parameters, existing for empirical purposes 
only. In reality, our knowledge of turbulence 
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in these stars is meager; the empirical curve-of- 
growth microturbulence is, in general, incon- 
veniently supersonic in some outer layers. 

The Mira characteristics (M. Querci, this 
volume), which will be discussed in the section 
Shock-Wave Gas Dynamics and Pulsation 
Theory, are more appropriate for a non- 
thermal concept. In this case, the velocities are 
no longer ad hoc parameters, and the shock- 
wave velocities that can exist in those stars are 
large enough to produce aerodynamic dissipa- 
tion. In any case, the complete nonthermal 
models for Mira and non-Mira-type giant and 
supergiant cool stars will be obtained in the 
future when a consistent solution of the 
dynamic-transfer problems will be found. 

Another important characteristic of the 
present models is that they are one-component 
models, and the probable existence of surface 
inhomogeneities is not taken into account. We 
know that large inhomogeneities exist in the at- 
mospheres of the Sun and of the late-type 
dwarfs. Giant cells are predicted to exist in the 
late-type giants (Schwarzschild, 1975), but to 
the present time, no satisfactory observational 
evidences of these cells have been advanced 
(see, for example, Johnson, this volume). One 
possibility for future development of multicom- 
ponent atmospheric models can be made by a 
simultaneous diagnosis of the hotter and cooler 
components of the atmosphere, for instance, 
by the Mg II emission and by the CO fun- 
damental vibration and rotation bands, respec- 
tively. (See Heasley et al., 1978, for an 
analogous study of Arcturus.) 

Two main regions are important for these 
stars: the photosphere and the chromosphere. 
In view of the observational evidence of the 
absence of a high-temperature plasma, no tran- 
sition regions and coronae are considered. The 
section Photospheres is devoted to the discus- 
sion of non-LTE effects in photospheres. In the 
section Chromospheres , we discuss the recent 
research on the chromospheres of these types 
of stars. In the section Shock- Wave Gas 
Dynamics and Pulsation Theory , we discuss 
elementary shock-wave and pulsation theories 


and their applications to Mira long-period 
variables. 

PHOTOSPHERES 

Our approach to the study of quasi-thermal 
photospheres will consist of relaxing the LTE 
hypothesis. In the subsection Kinetic 
Equilibrium Line Formation , we consider the 
non-LTE line-formation theory in photo- 
spheres. In the following subsection, we will 
discuss the metal ionization equilibrium and an 
application to a typical carbon-rich star photo- 
sphere. 

Kinetic Equilibrium Line Formation 

The main purpose of this subsection is to 
calculate a theoretical stellar spectra starting 
from a given atmospheric model representing 
the photosphere. This model consists essentially 
of a distribution in height of the kinetic tem- 
perature, gas and electronic pressures, and con- 
tinuum opacities. The validity of the Max- 
wellian velocity distribution of particles pro- 
voking the collisions is assumed. Non-LTE line 
formation using this assumption has received 
the name “kinetic equilibrium” (KE) line for- 
mation (Athay, 1972). 

Statistical Equilibrium. By definition, the 
statistical equilibrium equation relates the pop- 
ulation numbers of different levels of the atom 
(, n ) to the rates causing changes in this 
population. Note the important assumption 
that dn. /dt = 0. Considering an atomic 
model of N bound levels and a continuum (k), 
the statistical equilibrium equation can be writ- 
ten as: 



( 8 - 1 ) 
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i = i 

where P = R + C, and R and C are the 
radiative and collisional rates, respectively. The 
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upward rate is the excitation between two 
bound levels (/, /) or the ionization for a bound- 
free transition (/, k). The respective inverse pro- 
cesses are the deexcitation and recombination. 
The total number density of the atom is given 
by the particle conservation equation, 

N 

n tot = XI n ‘ + ”* ’ < 8 - 2 > 

1=1 


where n k is the number density of ionized 
atoms. Assuming that n (o{ is known and that 
the total rates can be calculated, the equations 
can be solved to find the populations n { and 

The radiative photoionization and recom- 
bination rates depend on the mean intensity of 
radiation J by the following expression: 


2hv 3 /c 2 

s = — - — 

L g j n i ' (8-5) 

*i n j 

g being the statistical weight. This last expres- 
sion can always be written in the form (using 
relation (8-1)): 
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$ v being the normalized absorption profile. 
This expression contains three fundamental 
physical processes: (1) the diffusion term, 
j 4 \Jjdv, (2) the source term or creation of 
line photons, eB', and (3) the sink term or 
destruction of line photons, e. 
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where a ( (v) is the photoionization cross sec- 
tion and v.. is the ionization threshold 

kl 

frequency. 

If J = B {T ), B being the Planck func- 
tion, the recombination rate can be calculated 
(see, for example, Mihalas, 1978) by 
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where n* / /n* k is given by the known Saha- 
Boltzmann equation. Knowing the collisional 
ionization rate, C /jt , a relation similar to (8-4) 
can be used to calculate the collisional recom- 
bination rate, C kr Radiative and collisional ex- 
citation rates can be calculated by classical rela- 
tions using the Einstein radiative coefficients 
and Maxwellian distribution of level popula- 
tions (see also Mihalas, 1978). 

For two levels / and y, the line-source func- 
tion is defined as: 


Transfer Equation. In reality, S L and J v are 
related by the equation of radiative transfer, 
measuring the change in distance of the inten- 
sity of radiation I ' , making an angle 0 (n = 
cos 6) with the normal to the surface: 

dl 

v =-MW • ( 8 - 7 ) 
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If we write the line-absorption coefficient as 
k L - k L (f> , where 0 = /<!> , we obtain the 

V O V V V o 

total monochromatic optical depth defined by: 
dr = (0 v + r o ) dr o , (8-8) 

where dr = k L dz and r = k c /k L . The 

o o o o 

transfer equation is then written as: 



(8-9) 


S being the total source function. 


S 

t 


0 S r + r S 

L o < 

+ r 


( 8 - 10 ) 


375 


where C refers to the continuum. The line- 
source function S L , as written in relation (8-6), 
admits the validity of the complete redistribu- 
tion of line photons (i.e., the emission profile 
is equal to the absorption profile). S L is then 
not frequency-dependent. The main goal of the 
KE line formation consists of solving 
simultaneously the statistical (8-6) and transfer 
(8-9) equations, specifying the boundary con- 
ditions and the geometry of the medium. 

Applications of KE Line Formation 
in LTE M-C Photospheres 

Very few works have been devoted to the KE 
line formation in photospheres of late-type 
giants and supergiants of the M-C type. This 
is not the case of the slightly hotter stars of type 
G and K, which have received much more at- 
tention. Regarding the cooler types, the main 
work is related to the abundance determination 
of lithium. The importance of this element as 
an indicator of the internal and atmospheric 
evolutionary state of these stars is well known. 
Due to the relative ease of creation and destruc- 
tion of this element, very large abundance dif- 
ferences can be found in these stars. These 
abundances vary from values smaller than the 
solar abundance to the extremely large values 
found in the “super lithium” C stars: WZ cas, 
WX Cyg, and T Ara. In the latter case, Li can 
no longer be considered as a trace element. 
However, abundances for these stars are not 
well known because they depend on solutions 
of the excitation/ionization problem described 
above, and these are not yet available. 

Lithium Line Formation of M-C Giant 
Photospheres. Recent work on Li abundance 
determination in M-type photospheres has been 
undertaken by Luck and Lambert (1982), us- 
ing improved observational data and based on 
model photospheres taken from the grid of 
Johnson et al. (1980). They used the complete 
linearization KE method code (Auer et al., 
1972) to calculate the non- LTE abundance of 
Li with an atomic model of four levels plus con- 


tinuum (Luck, 1977). Using an empirical solar 
line blanketing adapted to late-type giants, they 
found no non-LTE effects giving abundance 
differences for Li larger than 0.1 dex, sug- 
gesting that, because of the reduction of the 
ionizing field by line blanketing, non-LTE ef- 
fects are not important in the ionization balance 
of Li in early M giants and supergiants. This 
work is based on the resonance Li I line at 
X6708. It would be desirable for future Li 
studies to attempt to observe and use the subor- 
dinate Li I line at X6103. This weak line has 
been identified by A. M. Boesgaard (Merchant, 
1967) in early high-luminosity M stars. 

Concerning the carbon-enriched photo- 
spheres, de la Reza and Querci (1978) per- 
formed KE calculations for the neutral lithium 
lines. The method they employed for solving 
the radiative and transfer equations is the in- 
tegral flux divergence method (Athay and 
Skumanich, 1967; Athay, 1972). An atomic 
model of four levels and a continuum was also 
used. The photospheric LTE models were taken 
from the grid calculated by Querci (1972) and 
Querci and Querci (1975). The effective tem- 
peratures had values between 3800 and 3500 K, 
g = 10, C/O = 1.3. The photoionizing radia- 
tion field and the abundance were considered 
to be free parameters. The radiation field had 
the purpose of mimicking the effect of possible 
chromospheres and of simulating more or less 
blanketing. The main results were: 

1 . Li I is almost completely ionized in all 
models with F ff >3500 K. 

2. The ionization is increased in KE. 

3 . S L >B v (F) in the line-forming region, 
except for the case of the resonance 
X6708 line with an Li abundance equal 
to 10~ 8 . 

4 . Profiles are sensitive to Li abundance 
and the chromospheric radiation tem- 
peratures, but not to the possible 
presence of graphite grains in the 
temperature minimum region. 
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5 . Appreciable differences between KE 
and LTE were found for X6708 only, 
although the Li I line at X6103 is useful 
for distinguishing the effects of abun- 
dance from those of the chromospheric 
radiation field. 

6. In LTE, to produce lines of Li I as 
strong as those observed in the lithium- 
rich C stars without enhancing the Li 
abundance, would require T ff 2600 
K. 

Future Li abundance analysis should prob- 
ably take into account aspects that have not 
been considered in the past: heterogeneities and 
spherical symmetry. Recently, Giampapa 
(1984) found that an enhanced chromospheric 
activity, such as a solar plage, can diminish the 
equivalent widths (IV X ) up to a factor 2 due to 
the increased ionization. This agrees with the 
observed anticorrelation found between the in- 
crease of emission of Ca II K line and the W x 
of the X6707 Li line in the Sun. Special care 
must be taken when the Li I X6707 line is com- 
pared to lines of Ca or Na as, for example, to 
the neighbor line of Ca I at 6717 A. A proper 
KE analysis must also be made for these com- 
parison elements. 

Steenbock and Holweger (1984) investigated 
the ionization effect of the ultraviolet radiation 
field shortward of 3500 A, which controls the 
ionization balance of Li. The thermalization ef- 
fect of the H-Li atomic collisions was also 
estimated. In the case of a typical galactic disk 
late-type supergiant, the non-LTE effect 
measured by Alog e = log e - log e (LTE) can 
be equal to 0.57. This value is reduced by a fac- 
tor of 2 if line blocking is considered. Due to 
the low densities involved, the thermalization 
collisions with H atoms have a negligible effect. 
For late-type halo giants of smaller metallici- 
ty, the results are similar; however, due to the 
larger densities, the thermalization collisions 
can be somewhat more important. To examine 
the influence of W x on the Li abundance, a 
careful comparison must be made between the 
works of Steenbock and Holweger (1984) and 


Luck and Lambert (1982). Also, spherical sym- 
metry must be taken into account, particularly 
in the case of late-type supergiants. 

KE Line Formation of Molecules. The KE 
problem of line formation for molecules in 
stellar spectra can be significantly more com- 
plex than that for atoms. In fact, the analysis 
must contain all the rates of formation and 
destruction of molecules. In some cases, the 
destruction process, such as photodissociation 
of molecules, can be a minor process in the 
statistical equilibrium (i.e., the case, for exam- 
ple, of the important molecules TiO and SiO 
in M stars and CO and CN in C stars). These 
molecules have large dissociation energies, and 
the corresponding radiation for photodissocia- 
tion is located in the ultraviolet (X < 1500 A), 
which is very low in these stars, especially in 
the C stars. 

Other considerations that limit the ap- 
plicability of this analysis are the poor know- 
ledge of the cross sections involved and the 
numerical intractability of the calculation of a 
band formed by hundreds of molecular lines. 
An interesting and general work on these prob- 
lems is presented by Hinkle and Lambert 
(1975). Another work related to the formation 
of CO is that of Thompson (1973). The main 
conclusions of Hinkle and Lambert are as 
follows: 

1 . Because of the domination of the col- 
lisional process, LTE line formation 
(S L = B) is the appropriate form for 
vibration/rotation transitions occurring 
within the ground electronic state. 
(These are the transitions found in the 
infrared spectra.) 

2 . The equilibrium of the excited electronic 
state is radiatively controlled, especial- 
ly for late-type supergiants. These lines 
would be formed by the scattering pro- 
cess ( S L = \$Jjdv). For an optically 
thin absorption line, S L - J c v , where J c v 
is the local mean intensity in the con- 
tinuum. (These transitions are found in 
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the visible and near-infrared stellar 
spectral regions.) 

For the thin-line approximation, the equivalent 
widths will depend on the ratio, JJ B v , which 
can be very high in the outermost layers of 
stars. The equivalent widths will then be lower 
than those calculated in LTE. However, for the 
lines that are not optically thin, the line-source 
function depends on the scattering term, 
(<f> J dv, and the continuum source function. 
The line formation will be nonlocal in nature, 
and the complete KE equations must be solved. 
This complete KE calculation has been made 
by Mount et al. (1975) and Mount and Linsky 
(1975) for the violet CN band in the Sun and 
a Boo and by Carbon et al. (1976) for the fun- 
damental bands of CO in cool stars. All these 
investigations found significant departures 
from LTE. For late-type stars, S L < B has been 
found. Contrary to the thin case, these results 
give stronger absorption lines in KE than in 
LTE. 

The electronic bands, as is the case of CO, 
can be found in the ultraviolet. Ayres et al. 
(1981) showed evidence of the presence of 
fluorescent spectra of CO bands pumped by 
chromospheric emissions of O I at X1380 and 
C I at XI 600. These bands appear as weak emis- 
sions near the C II (XI 335) and C IV (XI 550) 
features, which are normally indicators of the 
presence of a hot gas at 10 5 K. These authors 
then call for caution in the interpretation of 
these features in red giants. 

Ionization Equilibrium in Photospheres 

Non-LTE effects have always been sus- 
pected to exist in the atmospheres of late-type 
giant and supergiant stars due to their low total 
densities. However, rather than the total den- 
sity, it is the free electron density which is im- 
portant, because of its contribution to the ex- 
citation/ionization collisional process and to 
the continuous opacity. The goal of the study 
of ionization equilibrium is to determine the 
real distribution and effects of the electron den- 
sity produced by the electron donors. In cool 


stellar photospheres, instead of H, which is the 
fundamental electron donor in chromospheres, 
the most abundant metals with low-ionization 
potentials are the principal sources of electrons. 
These are principally K, Na, Al, Mg, and Ca. 

In this respect, very little study has been 
directed to these interesting problems in the 
photospheres. We can cite only two investiga- 
tions: the theoretical work of Auman and 
Woodrow (1975) and the empirical work of 
Ramsey (1977, 1981). Ramsey compared the 
observed ratios of the equivalent widths of the 
[Ca II] X7323 and Ca I X6573 lines with those 
calculated by LTE model atmospheres of late- 
type supergiants. A certain discrepancy of these 
ratios begins to be detected for models cooler 
than 4250 K. This could be interpreted as 
evidence of non-LTE effects. However, as Luck 
and Lambert (1982) show, the continuum 
placement in these spectral regions is difficult, 
suggesting an inadequacy of the apparent local 
continuum of Ramsey. Nevertheless, we must 
not be far from reality, and as the theoretical 
work of Auman and Woodrow show, there are 
clear indications of non-LTE effects in high- 
luminosity cool stars. These authors iterative- 
ly solve the atmospheric equations with only the 
statistical equilibrium of the metallic electronic 
donors mentioned above. The simplifications 
are: (1) no radiative transfer is assumed, (2) the 
bound-bound transitions are treated in radiative 
detailed balance, and (3) no atomic or molec- 
ular (with the exception of H 2 0) blanketing is 
taken into account. Even with these simplifica- 
tions, their work was a fundamental step in this 
type of research. 

To discuss the basic features of the ioniza- 
tion equilibrium, we consider the simplest case 
of a one-level atom with a continuum. The sta- 
tistical equilibrium (see Equation (8-1)) is then: 

",<*/* + C ,P = n k (R kl + C„). (8-11) 

Using Equations (8-3) and (8-4) and the 
Saha-Boltzmann relation, Equation (8-1 1) can 
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be solved to give the value of the departure 
coefficient, b y which is given by: 



When the atoms are predominately neutral, 
almost all atoms are at level 1; therefore: 

b { ~ n\/n k , (8-13) 

and when atoms are strongly ionized, n k ~ 
n* k , b is reduced to: 

b { ~ /ij/zi* . (8-14) 


Actually, the mean intensity, 7 v , appearing 
in Equation (8-11) is also a function of b v and 
the expression relating J y to b can be found by 
simultaneously solving the statistical and 
transfer equations (KE). Since the line- 
absorption coefficient is also a function of n ( 
and b v an iterative process must be used to 
solve the equations for the values of 
b x . Due to the low electron densities existing 
in the atmospheres of late-type giants, the value 
of b~ l depends exclusively on the ratio J v /B y . 
(For an explicit form of Equation (8-12), see 
Avrett and Loeser, 1969.) This ratio may be 
very large in the outermost regions due to the 
large values of the Rayleigh scattering coeffi- 
cient compared to the other continuous opac- 
ities, especially in the ultraviolet. In the deeper 
layers, the ratio, J v /B yS converges to unity, 
giving a departure coefficient b { equal to 
unity. 

A practical form to express the ionization 
equilibrium is (Auman and Woodrow, 1975): 



where N + , TV, and N are the ion, neutral, and 
electron densities. U and U* are the partition 
functions assuming KE and LTE, respectively. 
It is easy to see that, in general, b~ x ~ U*/U. 


Assuming the neutrality of the plasma, the new 
electrons can be calculated by means of the new 
ion density. This “new” model can be iterated 
to obtain new values of J v which can be used 
again in the KE calculation. This iterative pro- 
cess can be used until convergence in the depar- 
ture coefficients is obtained. 

To illustrate the behavior of the b { values 
in model photospheres, we have calculated the 
KE for neutral atoms of K, Na, and Mg (for 
5, 4, and 5 bound levels, respectively). These 
calculations consist of finding the first values 
of U*/U , using the mean intensities J v from 
an LTE molecular-blanketed model photo- 
sphere. The model chosen is that of Querci 
(1972) and Querci and Querci (1975) corres- 
ponding to T eff /g/ratio of C/O equal to 
3400/0.1/3.26. Because no atomic line or grain 
blanketing was considered in those models, the 
J v values, especially in the ultraviolet, are 
maximum values. The U/U* obtained, then, 
expresses maximum departures from LTE. It 
is important to mention, however, that subse- 
quent iterations with model calculations will 
reduce, even more, the b } values. These results 
are presented in Figure 8-1. 

We have experimented numerically to gage 
the sensitivity of U*/U to the changes of J y 
values at different continua. We find that a 
reduction in J by a certain factor for the first 
excited levels of all atoms produces a reduction 
of approximately the same factor in the values 
of U*/U. These results are also sensitive to the 
second excited state of Mg I (as also found by 
Auman and Woodrow, 1975). We have ex- 
amined the abundance effects by arbitrarily 
changing the model abundances for K and Na. 
Although there are no significant changes in the 
J v values of the LTE model, the KE popula- 
tions of the bound levels change, giving larger 
departures for smaller abundances. 

The converged values of Auman and 
Woodrow do not show significant temperature 
changes from the initial models. The most im- 
portant changes are found in the electron den- 
sity for the models: T ff = 2500 K, log g = 
0.0, and T ff = 2000 K, log g = -1.0, in which 
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Figure 8-1. Behavior of the ratio of the initial 
partition functions with depth in the photo- 
spheres of a carbon-star model and for the prin- 
cipal electron donors , K, Na, and Mg. m is the 
column mass. A scale of depth for the reference 
wavelength, \8063, is also shown. 

the electron pressures increase in the outer 
layers, compared to LTE models, by factors of 
14 and 60, respectively. Only slight changes are 
found in the continuous opacity (in general, less 
than 5 percent). This is due to the fact that, in 
the outermost layers where the non-LTE effects 
are larger, the continuum opacity is controlled 
by Rayleigh scattering by atomic and molecular 
hydrogen (especially in the ultraviolet), which 
is independent of Ne. Of course, H~ is an im- 
portant source of the continuous spectra in 
these stars, and it depends on Ne. However, 
H~ appears to be important only in deeper 
layers in which the non-LTE effects are smaller. 
It is also important in this respect that a com- 
plete calculation of the H" equilibria be per- 
formed following the suggestions of Lites and 
Mihalas (1984). 

We conclude that non-LTE effects found 
for the carbon-star astmospheres shown in 
Figure 8-1 are large. These maximum values, 
however, are partial and will be reduced by 
iteration with the model. We hope that future 
models will contain blanketing by atomic lines 


(and grains?) to enable construction of reliable 
non-LTE blanketed atmospheres. 

CHROMOSPHERES 

The chromospheres of noncoronal stars 
such as M and C high-luminous stars are 
regions that lie immediately above the 
photospheres and in which a large part of the 
nonradiative energy available heats and ex- 
pands the gas in order to produce an extended 
warm region. It is in these regions that the mass 
loss begins to take place. 

Nevertheless, it is precisely in this type of 
star that real atmospheres seem to be much 
more complicated than in our simple models. 
That is due to the presence of a number of 
physical processes involved (e.g., abundance 
and grain effects, weak and strong shock 
waves, magnetohydrodynamic waves, turbu- 
lence, and possible large convection effects). 
Although a realistic energetic budget will be dif- 
ficult to calculate under these conditions; it is 
the only way to understand why some stars 
have different kinds of chromospheres. (See 
Linsky, 1984a, for a recent review on 
chromospheres of non-Mira M and hotter G- 
and K-type stars.) 

Hydrogen Ionization 

In photospheres, ionization of abundant 
low-ionization metals was important, whereas 
in chromospheres, the ionization role is played 
by hydrogen. In fact, even at the modest tem- 
peratures, between 4000 and 10000 K, which ex- 
ist in these zones, hydrogen will be ionized. 
Ionization of helium is less important, since it 
will be produced only at higher temperatures. 

The construction of chromospheres must 
then calculate simultaneously the H line and 
continuum formation, together with the at- 
mospheric equations using methods, for exam- 
ple, like those developed by Avrett (for a re- 
cent paper, see Avrett, 1985). H continua can 
also be treated by a method developed by 
O’Brien (1980) using the code (LINEAR) of 
Auer et al. (1972; see also Auer, 1973). This last 
procedure can be modified to treat H lines also. 
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Regarding H, another interesting problem 
that should be mentioned is the puzzle of the 
H 2 quadrupole line in M and S stars (Tsuji, 
1983) and in C stars (Goorvitch et al., 1980; 
Johnson et aL, 1983). The strongest infrared 
H 2 quadrupole line is not detected in non-Mira 
M stars, contrary to the theoretical expectation 
which predicts an increasing strength of this line 
as we go to cooler stars (Tsuji, 1983). This H 2 
feature is clearly observed in S-type Mira stars 
(Hall and Ridgway, 1977; Tsuji, 1983). In C 
stars, the line is present in Mira variables. 

Tsuji (1983) and Johnson et al. (1983) 
describe the principal physical scenarios that 
could explain the presence or absence of this 
line. Even if the problem between M and C 
stars appears to be similar, it is not clear 
whether the explanation will be the same for 
both types of stars. 

Lower surface temperatures produce a 
stronger H 2 line if we consider the thermal 
sensitivity of the line. Spherical symmetry can 
produce these lower temperatures (Schmid- 
Burgk et al., 1981). This effect could thus ex- 
plain the presence of H 2 in S stars. Tsuji 
(1983), however, favors the explanation of 
lower surface temperatures in S stars in respect 
to M stars (differences between 150 and 350 K), 
produced by the more efficient cooling effect 
of CO in S stars. Nevertheless, the problem 
concerning the non-Mira M stars remains open. 
Unless a general revision of the temperature 
scale is made, other mechanisms, such as the 
presence of a chromosphere, particularly on 
later types, could be invoked to explain the 
destruction of the H 2 molecules. 

It is possible that the puzzle has an explana- 
tion other than thermal sensitivity for the case 
of C stars. Johnson (see Chapter 7) suggests the 
possibility of hydrogen deficiency in C stars. 
This suggestion is strengthened by an impor- 
tant characteristic of C stars (i.e., the absence 
of Balmer lines in N carbon stars (Yamashita, 
1972, 1975)). Another fact that suggests an H 
deficiency is the weak observed emission peak 
of H" at 1.6 fx m in N stars as compared to that 
predicted by models (Goebel and Johnson, 


1984). These authors propose that, if H defi- 
ciency is real in N carbon stars, this can be ex- 
plained by an episodical ejection of an H-rich 
shell. That conclusion has important conse- 
quences for the evolutionary study of these ob- 
jects and the chemical evolution of galaxies. It 
is important in this respect, however, to cal- 
culate the H“ equilibrium in the most complete 
way, following the treatment given by Lites and 
Mihalas (1984). Perhaps we could explain the 
absence of Balmer lines in N stars by the 
mechanism put forward by them, in which 
charge neutralization processes decrease excited 
H-level populations in favor of protons and 
H' ions affecting the H ionization and the 
Balmer lines. 

Bifurcated Chromospheres 

The behavior of such important heating and 
cooling sources as H“, CO bands, and metallic 
lines in stellar gases at the minimum solar 
temperature region, r min ~4000 K, has in- 
troduced new ideas on the structure of the solar 
chromosphere (Ayres, 1981). For the latest 
discussions, see Ayres (1985) and Avrett (1985). 
We think it will be interesting to check whether 
these concepts can be applied to future studies 
of detailed late-type star atmospheres. 

First, let us examine the mechanisms: from 
one side, CO cooling is important for temper- 
atures lower than 4000 K. For T > 4000 K, the 
molecule dissociates, and its cooling effect 
diminishes. (The role of CO as a cooling source 
in late-type stars has been discussed initially by 
Johnson, 1973.) On the other hand, at this 
r min ~ 4000 K region, H“ is contrary a heating 
source due to photodetachment. H’ will be 
converted in a coolant only at higher tem- 
peratures near 5000 K. Bearing this in mind, 
we can now describe the principal physical idea 
(Ayres, 1985). For temperatures near 4000 K 
in a stationary atmosphere, the cooling proper- 
ties of CO can act against the heating produced 
by H' and the deposition of nonradiative 
energy, maintaining a cool gas at T ^ 4000 K 
up to higher zones which were classically de- 
fined as belonging to the chromosphere. 
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However, in some supposedly restricted areas 
of the photosphere surface (e.g., magnetic flux 
tubes), the deposition of nonradiative energy 
can be strong enough to destroy the CO 
molecule and produce a temperature inversion 
that forms a localized chromosphere, for which 
the classical 7 . concept can be applied. This 
stationary atmosphere will then be a two- 
component thermal atmosphere in which two 
zones — one hot, the other cold — will coexist at 
the same height. This type of bifurcated model 
can also be explained by nonstationary at- 
mospheres (Muchmore and Ulmschneider, 
1985) in which waves produced hot compress- 
ed zones and rarefied cool zones cooled further 
by CO radiation. Kneer (1983) has also examin- 
ed the possibility that the thermal instability of 
optically thin infrared vibration/rotation bands 
of CO could produce chromospheres of late- 
type stars, proposing an interpretation of the 
known Wilson-Bappu effect (Wilson and Bap- 
pu, 1957). 

One of the purposes of these bifurcated 
models is to reproduce the ionic emission as Ca 
II K simultaneously with the radiation of the 
CO infrared band. The empirical basis of these 
models is presented in Ayres (1985). Future 
development will depend, as Ayres explains, on 
a detailed treatment of non-LTE effects of CO, 
especially for ultraviolet electronic transitions. 
Metallic blanketing will be important because 
non-LTE effects will not only decrease the cool- 
ing effects of metals by a large factor but will 
also increase their effectiveness in depth by a 
similar factor (Kalkofen, 1985). 

We would like to emphasize that it will be 
interesting to investigate the application of 
bifurcated models to late-type stars for which 
simultaneous observations of metallic emission 
lines and CO bands have been made. 

Chromospheres of M-Type Stars 

Observational Evidences. The research on 
chromospheres in M-type giants and 
supergiants has been based on the presence of 
some emission lines (see, for example, M. Quer- 
ci, this volume) like the Mg II and Ca II reso- 


nant lines (Basri et al., 1981; Hagen et al., 1983) 
and Fe II lines (Boesgaard and Magnan, 1975; 
Boesgaard, 1979). Other lines also appearing 
in emission are the hydrogen Balmer lines (Jen- 
nings and Dyck, 1972; Hagen et al., 1983) and 
the He line (Wilson, 1957, 1982). However, the 
hydrogen lines can be explained by other 
mechanisms (e.g., shock waves). In the 
ultraviolet, research is also concentrated on the 
C II lines. 

Dyck and Johnson (1969) and Jennings and 
Dyck (1972) first showed an anticorrelation be- 
tween the H and K lines of Ca II and the infra- 
red excess and polarization, characterizing the 
presence of dust grains in their atmospheres. 
This anticorrelation was interpreted by Jen- 
nings (1973) as due to a quenching effect by 
dust grains. These grains would radiate the non- 
thermal energy, preventing the formation of a 
chromosphere. This anticorrelation has also 
been found by Hagen et al. (1983) in a relatively 
large sample of M stars. However, the results 
of Hagen et al. show a difference in the sense 
that the real anticorrelation of Ca II emission 
is not with circumstellar dust, but appears to 
be with a dust-to-gas ratio (the degree of grain 
condensation). The core emission of the Ca II 
lines disappears in the coolest stars having 
strong dust-to-gas ratios. That is the case, for 
example, for the stars, HD 207076 (M7 III) and 
R Dor (M8 III). We can infer that these stars 
probably do not possess a chromosphere. How- 
ever, any conclusion must be analyzed in detail 
because R Dor also shows Mg II in emission, 
which is normally considered as a chromo- 
spheric indication. Hagen et al. suggest a possi- 
ble shock origin for this Mg line. This shock 
origin is also suggested for the Balmer emission 
lines that appear only in stars with high dust- 
to-gas ratios (Jennings and Dyck, 1972; Hagen 
et al., 1983). 

Extension of the Chromospheres. Recent work 
based on the ultraviolet (UV) C II emission lines 
has permitted an empirical estimate of the 
temperatures, electronic densities, and 
geometric extensions of chromospheres of 
luminous stars of spectral types between G9 and 
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M3. The C II lines involved are the resonance 
triplet (UV 1) at 1335 A and the multiplet (UV 
0.01) at 2325 A. Stencel et al. (1981) showed 
that, for electronic densities between 10 7 < N 

e 

< 10 9 cm' 1 , the line ratios of the multiplet, 
UV 0.01, can be used as an indicator of N . 
The electronic temperatures are obtained by 
comparison of the observed ratios of the UV 

I and UV 0.01 total fluxes, together with a 
theoretical calibration of this ratio with N and 

e 

T e like that of Hayes and Nussbaumer (1984). 
Brown and Carpenter (1985) introduced this 
technique to obtain the temperatures of the C 

II line-emitting regions in several luminous late- 
type stars, assuming that the radiation arising 
from both multiplets came from a region of 
constant N g . For M stars, the chromospheric 
values are N - 10 8 cm" 3 and 7000 K£. T < 

e e — 

9000 K. The most interesting estimations are, 
however, those related to the sizes of these 
chromospheric C II emission regions. The most 
recent results are those of Carpenter et al. 
(1985), obtained by using the observed total 
fluxes of multiplet UV 0.01, the above- 
mentioned values of N , T 9 the photospheric 
radii, R*, and the angular diameter, <£*. The 
mean thickness of the emitting region, L, is 
calculated by these authors by the relation: 


/ 23 25 

L = 

A E 


(8-16) 


where A is the solid angle subtended by the 
chromosphere and E is its volume emissivity 
given by: 


£ = 7^cn c 12 


(8-17) 


The number density of C II can be calculated 
by: 
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— N. , (8-18) 


for which the expression of fractional ioniza- 
tion of C and H and C abundance must be 


known. The solid angle is expressed in terms 
of the chromospheric radius, R , and the values 
of R* and </>*. The principal unknowns, L and 
R , are then related, and an iterative technique 
is used. These estimations, although subject to 
uncertainties, as values of N and the C abun- 
dance, clearly show that late-type giants with 
coronal regions (G9 and K0) have very thin 
chromospheres, with R < 0.1 percent of the 
photosphere. Noncoronal stars like M giants 
and supergiants have instead, on the average, 
chromospheric emission regions extending out 
to 2.5 photospheric radii. In the case of the M 
star, a Ori, for example, this method gives a 
chromospheric radius smaller by a factor of 2 
than that found by Newell and Hjellming (1982) 
using radio observations. Of course, one must 
consider that the assumption of constant N g is 
quite simplistic. 

Thus it appears that these red high- 
luminosity stars contain a vast chromosphere 
extending to several stellar radii many times 
larger than their isothermal pressure scale 
heights. Mass loss is bound to be happening in 
these regions. Kinematically speaking, the 
chromospheres could, in a way, be physically 
detached from the photospheres (Goldberg, 
1979). How can these extended chromospheres 
be produced? Certainly, the mechanism of ther- 
mal evaporation, which works in a hot stellar 
wind, will be inoperative here, due to the low 
temperatures and thermal pressures involved. 
In addition, it is not clear whether radiation 
pressure impinging on molecules (Maciel, 1976, 
1977) or La radiation (Wilson, 1960; Haisch et 
al., 1980) could produce those extended regions 
by gas levitation and expansion. This will re- 
quire further study. Other pressures of tur- 
bulent or magnetic origin can also be invoked. 
Stencel (1982b) mentioned that an rms tur- 
bulent velocity of the order of 50 km s" 1 would 
be necessary to support material with a densi- 
ty of 10 8 cm" 3 at 5 stellar radii. A magnetic 
pressure of a modest magnetic surface field 
could produce this support (see Goldberg, this 
volume). 

Line Formation in Extended Atmospheres. If 
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chromospheres extended to several photo- 
spheric radii, spherical geometry must be taken 
into account in any realistic treatment of line 
formation. In addition, chromosphere expan- 
sion effects must be considered if mass loss is 
taking place somewhere in these regions. 
Because wind outflow velocities of the order of 
10 to 50 km s' 1 are only a factor of 2 larger 
than typical microturbulence values in late-type 
giant stars, the Sobolev approximation that is 
valid for strong winds will probably be of no 
value. Recent numerical applications solve the 
transfer equation in a comoving frame for a 
spherical atmosphere in expansion, assuming 
a partial redistribution (PRD) of line photons 
(Drake and Linsky, 1983; Drake, 1985 (dif- 
ferential equation method); Avrett and Loeser, 
1983 (integral equation method)). 

For the photospheres, Equation (8-6) admit- 
ted the validity of the complete redistribution 
(CRD) function, in which the emission-line pro- 
file, is equal to the absorption profile, 4^. 
This is no longer the case in the interpretation 
of strong resonance chromospheric lines. Here 
the PRD approximation represents an in- 
termediate situation between the coherent scat- 
tering and the CRD (in which, contrary to the 
coherent case, no frequency relation exists be- 
tween the incoming and outcoming photons). 
A recent numerical solution is presented by 
Hubeny and Heinzel (1984) on the effects of 
PRD on subordinated lines. 

The total actual probability of a photon of 
frequency v being absorbed and reemitted at v ' 
is: 


R(v,v') = y R u (v,v’) 

(8-19) 

+ (1-7) (*0 4 > (O , 


where R u is the angle averaged probability, 
and 7 is the coherence mixing fraction given by: 


7 = 


r rad^I 
r rod + r i +r E 


( 8 - 20 ) 


where R is the angle averaged probability, 
and 7 is the coherence mixing fraction given by: 


inelastic deexcitation collisions, and elastic col- 
lisions (i.e., Stark and Van der Waals col- 
lisions). 

Noncoherent scattering (CRD) is valid in the 
Doppler core, whereas coherent scattering is 
important in the inner wings. Absorption by 
continuum processes is important in the ex- 
treme wings. The difference between PRD and 
CRD has been discussed by Basri (1980) in 
static atmospheres and by Drake and Linsky 
(1983) in expanding chromospheres. The results 
of Drake and Linsky show, for example, that 
the increase of flow velocity and geometrical 
extension gradually transforms a double peak 
Mg II line into an approximate P Cygni Mg II 
line. Drake (1985) has described and applied the 
PRD techniques in a semiempirical way to in- 
terpret the Mg II line in Arcturus. However, 
no unique atmospheric model that explains the 
line profile has been found. As Drake men- 
tions, it is important to make future tests of 
cooler giants that show other chromospheric 
line indicators. Also in this respect, it is instruc- 
tive to compare the results obtained with those 
produced by pure theoretical models. For an 
extensive review paper on PRD stellar applica- 
tions, see Linsky (1984b). 

Model Diagnosis. Two different approaches are 
generally made to model chromospheric 
temperature, density, and velocity distribu- 
tions. One is to calculate the KE of important 
resonance lines, usually Mg II and Ca II for 
various assumed chromospheric models, and 
then select the model which best fits the ob- 
served profiles. This semiempirical method has 
been applied by several workers in F-, G-, and 
K-type stars (see, for example, Linsky, 1980, 
1981). Some elements of this methodology are 
presented in the section Chromospheres of Car- 
bon Stars. The other approach is a theoretical 
one that consists of predicting a chromospheric 
model from an adopted nonradiative heating 
mechanism. This theoretical treatment has also 
been applied in the study of F-, G-, and K-type 
stars (Ulmschneider, 1979; Musielak and Sikor- 
ski, 1981). An important general theoretical 
work on this field is that of Athay (1981). 
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Very few detailed chromospheric studies of 
M-type stars exist, and to our knowledge, only 
the supergiant a Orionis (M2 lab) has received 
special attention. Lambert and Snell (1975) first 
examined whether a simple nonextended 
constant-temperature chromosphere could con- 
tribute to the infrared excess observed in a Ori 
and in the semiregular variable, W Hydrae 
(M7-9). This contribution would be important 
only at wavelengths for which the silicate dust 
is a poor emitter. Isothermal chromospheres of 
5000 and 8000 K were proposed for a Ori and 
W Hya, respectively. However, this kind of ap- 
proach requires a detailed knowledge of the 
sometimes badly known dust properties. In ad- 
dition, those chromospheres must now be 
placed in a huge CS region in which grains ex- 
ist. New problems are now emerging as to 
determining the real distribution of grains in 
space and time in those huge atmospheres. For- 
ming dust grains near the surface of a Ori also 
appears to be difficult (Draine, 1981). 

Another approach in the study of 
chromospheres of M stars has been to analyze 
the velocity fields in a Ori (Boesgaard and 
Magnan, 1975; Boesgaard, 1979) and in the 
M2. 5 giant, 0 Pegasi (Boesgaard, 1981); these 
studies were devoted to the formation of their 
strong Fe II emission lines in the 3100 A region. 
The transfer problems are solved in a spherical 
moving atmosphere covering an extended 
region around the photosphere of the star. The 
line-source function is not specifically solved, 
but rather, it is parameterized in order to repre- 
sent a chromospheric rise of the temperature. 
Finally, the optical depths are adjusted, by 
varying the velocity fields, to fit the observed 
profiles. One important result for a Ori is that 
the Fe II lines in this model are formed in a 
region that shows an infall of matter. The 
velocity changes from 15 km s" 1 at 1.8 R* (i?* 
is the stellar radii) to 60 km s" 1 at the 
photosphere. In the giant, 0 Peg, the Fe II emis- 
sion lines show a rather uniform outflow begin- 
ning at zero velocity at the photospheric sur- 
face and attaining 25 km s" 1 at 2 i?*. 

A more direct semiempirical approach to the 
chromosphere of a Ori has been undertaken by 


Basri et al. (1981), who fit KE partial redistribu- 
tion calculations with the observed profiles of 
Mg II and Ca II resonant lines. The computed 
models are nonextended chromospheres in 
hydrostatic equilibrium, resulting in a plane- 
parallel atmosphere. These authors could not 
find a homogeneous chromosphere for a Ori 
that fitted both the Ca II and Mg II profiles. 
This preliminary model is shown in Figure 8-2. 

Zarro (1984) has recently studied the 
possibility that the profiles of Balmer lines, 
such as Ha, could be chromospheric indicators 
in late-type giants. However, as we will see later 
in the section Shock - Wave Gas Dynamics and 
Pulsation Theory , H-lines can be formed in 
some M giants by processes other than 
chromospheric ones (i.e., strong shocks). Also 
in this respect, in N irregular carbon-type stars, 
the Balmer lines are extremely weak 
(Yamashita, 1972, 1975). 

A different approach (theoretical) has been 
taken by Hartmann and Avrett (1983), who 
construct a model for a Ori by including the 
Alfvenic wave pressure to produce an extended 
chromosphere. The momentum and energy 
equations of the wind are solved for a radial 
steady spherical symmetric flow. The 
temperature, density, and velocity of the flow 
is then later computed to attempt to reproduce 
some observed features in a Ori. The tem- 
peratures are those resulting from the energetic 
balance of wave heating, adiabatic expansion, 
and radiative cooling. 

The principal difficulty of this mechanism 
stems from our limited understanding of the 
dissipation of these waves. The geometrical 
scale in which this dissipation occurs is 
represented by the damping length, L. Con- 
sidering L to be constant, the physical situation 
is very sensitive to the values of X = L/i?*. If 
X > 1 , the theoretical terminal velocity of the 
wind greatly exceeds the observed values, 
whereas for X < 0.85, no mass loss occurs, 
Hartmann and Avrett (1983) choose an in- 
termediate constant value X = 0.9 for their 
calculations. The two computed models are 
shown in Figure 8-2. The 4 ‘high-expansion' * 
curve corresponds to a model giving a mass loss 
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Figure 8-2. Temperature versus height chromospheric models of the supergiant M star, a Orionis. 
The nonextended semiempirical atmosphere ofBasri et al. (1981) is shown, together with the theoretical 
extended models by Alfven pressure waves (high and low expansion) of Hartmann and Avrett (1983). 
The photospheric zone is the same for all models. 


equal to 1.4 x 10' 6 A/© yr' 1 . The expansion 
velocities (V ) are of the order of 20 km s' 1 
and are larger than the turbulent velocities 
(V A. The other model, “low-expansion” 

v turb ' 

with V < V . , corresponds to a static 

exp turb 

energy equation. The photospheric model is the 
same as that of Basri et al. (1981). Both the 
high- and low-expansion models produce 
similar line and radio continuum emission in 
reasonable agreement with observations. 
Discrepancies with observations appear only 
when observed and computed profiles are com- 
pared. Theoretical profiles are much narrower 
than those observed. Predicted asymmetries are 
also too large. (The low-expansion model, 
however, produces nearly symmetric profiles.) 
In conclusion, larger V turh and smaller V gxp by 
factors of 2 and 3, respectively, are necessary 
to match the observed data, suggesting a 


“quasi-static” extended chromosphere for 
a Ori. 

Holzer et al. (1983) have criticized the use 
of a constant damping length. They believe that 
a more realistic interpretation will be one in 
which L is nonlinearly related to the state of 
the atmosphere (temperature, density, ioniza- 
tion, and radiative losses). They show that, if 
an ion-frictional damping is used in which the 
wave dissipation is produced by charge-transfer 
collisions between ions and atoms, it will give 
a thermal structure very different from that ob- 
tained by assuming constant L. These new 
temperatures are also the result of balance be- 
tween heating by dissipation and cooling by 
radiation obtained in a self-consistent way. 
Holzer et al. (1983) also raise objections to the 
Alfven wave mechanism as a general way to 


386 


produce cool massive low-speed winds for low- 
gravity stars. For them, any efficient 
mechanism (Alfven waves or something else) 
must add the major part of the available energy 
in the subsonic region, without affecting the 
asymptotic flow speed u^. ( u ^ is much less 
than the escape velocity in low-gravity stars.) 
A small remaining energy will then be neces- 
sarily added to the flow in the supersonic region 
to guarantee the real escape of the wind. 

Future modelers of a Ori must consider that 
time-dependent phenomena can be important 
input data. Although a Ori is a long-period 
variable ( P - 21 10 days), some rapid variations 
on a monthly scale have been observed 
(Goldberg, 1979). Other time-dependent 
phenomena could also be associated with possi- 
ble large convective cells (Schwarzschild, 1975) 
or other mechanism, producing sporadic ejec- 
tion of matter with velocities generally inferior 
to the escape velocity, so that infalling gas will 
be common, supporting perhaps the infall ac- 
celeration suggested by Boesgaard (1979). In 
addition, we cannot avoid the possible 
coexistence of hot and cool material produc- 
ing an inhomogeneous multicomponent model. 

Chromospheres of Carbon Stars 

Although the presence of some Fe II emis- 
sion lines in the violet appear to indicate the 
existence of chromospheres in N-type carbon 
stars, confirmation in certain stars (N stars) has 
come only recently through IUE observations 
of the 2200 to 3200 A interval (Johnson and 
O’Brien, 1983; Querci et al.,1982; Querci and 
Querci, 1985), and in R stars by Eaton et al. 
(1985). 

The stars observed by Johnson and O’Brien 
(1983) were all warm N-type stars, and none of 
them showed a strong violet depression. In fact, 
their ultraviolet spectra, beginning at 3200 A, 
show a continuation of the photospheric spec- 
trum with several absorption features down to 
wavelength X2850, where the spectra appear to 
change and present emission lines. Here, two 
clear emission features are present that corres- 
pond to the Mg II and C II lines at X2800 and 


X2325, respectively. Because of the low resolu- 
tion of IUE, the Mg II h and k lines are not 
resolved. The ratios /(Mg II Vf hol obtained are, 
on the average, 100 and 25 times smaller than 
the respective ratios of K and early M stars. 
However, the ratio /(C \\)/f bo/ is only 2 to 4 
times smaller than that of M giants. This dif- 
ference is partially a consequence of the 
enhanced carbon abundance in carbon stars 
relative to M-type stars. The ratio /(Mg II)//(C 
II) could be considered, in a way, to be indepen- 
dent of any superposed absorption. Even in this 
case, the latter ratio is 4 to 5 times smaller than 
that of M stars. From these results, Johnson 
and O’Brien suggest that chromospheres, 
weaker than those of M stars, must exist in 
these observed N-type stars. 

How can the absence of light detection from 
other target stars be explained? Some of them 
have been integrated with exposure time up to 
5 hours. Querci et al. (1982) explain this in two 
ways: one is the possible absence of a chromo- 
sphere. (This chromosphere could also be vari- 
able in time.) The other is the existence of a sup- 
plementary “ultraviolet depression.” It is not 
clear if this supplementary opacity is or is not 
related to the classical violet depression that ex- 
ists in some C stars and which begins at - 4200 
A. At this point, it is important to note that 
the detection of light depends on the brightness 
of the star at U. Only the light of the brightest 
U stars like TX Psc and TW Hor were detected, 
whereas faint U stars such as YCVn and WZ 
Cas were not detected. This simple fact in- 
dicates that this is due to ultraviolet opacity of 
molecular and/or grain origin. (For a discus- 
sion of grains, see Lefevre, this volume.) 

Simple atmospheres can be formed by a 
photosphere and empirical plane-parallel ad- 
hoc chromosphere. When the photosphere is 
known, the chromosphere can be calculated us- 
ing the following relations: 

1. The hydrostatic equilibrium expressed 
by: 

P-P 0 = (8-21) 
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where P and m are the total gas 

o o 

pressure and the column mass at the top 
of the chromosphere. 

2 . The total pressure given by the sum of 
partial pressures of the most abundant 
elements (in number N g[ ) 

N el 

P = n H ^ J (A { +N e )KT , (8-22) 

/ 

where is the abundance of the ele- 
ment /. 


3 . The charge neutrality, 



/ 


where N is the proton number, and 
each first stage of ionization is given by 
the function, f r between neutral 
elements and ions of the element, /. For 
H , the statistical equilibrium must be 
solved to obtain the departure co- 
efficient, 6 1# for the ground level. For 
H and other metals, ionization can be 

e 

given by the Saha law. 

We have four unknowns: N H , P, N g , and 
TV for Equations (8-21), (8-22), and (8-23). 
The fourth equation must necessarily result 
from the energy-momentum conservation equa- 
tion at any depth of the chromosphere. One 
way to replace it is to obtain an arbitrary rela- 
tion between temperature and mass. One sim- 
ple way is given by: 

T = a log m + j3 , (8-24) 

where a and /3 are constant coefficients deter- 
mined from the boundary conditions: T q at 
the top of the chromosphere and T mjn at the 
minimum temperature region corresponding to 
the top of the photospheric model. All these 


equations can be solved by a rapid iterative pro- 
cess to produce the chromosphere. 

The energy boundary conditions must be es- 
tablished for the fundamental variables, P and 
T: P , T and P and T . The choice of a 

o o nun nun 

smaller value of P q will produce a more ex- 
tended chromosphere for a given temperature 
gradient. 

Using this method and keeping in mind the 
exploration of possible chromospheres of car- 
bon stars, we have calculated some linear 
chromospheric extensions to the photospheric 
models of Johnson (1982). For these calcula- 
tions, we have chosen two different carbon- 
enriched model photospheres with very dif- 
ferent gravities and different effective tem- 
peratures. The main parameters (T, ff /log 
g/ratio of C/O) are 3500/0.0/1.05 and 
3000/2.0/1.02. For a chosen temperature dis- 
tribution, the model will contain the calculated 
values of the column mass, the electron den- 
sities, and the departure coefficient, b ( , of the 
hydrogen. In the photosphere, b { will be equal 
to unity. We have also calculated an ad-hoc 
depth-dependent microturbulence (increasing 
outward) by means of the following relation: 

£ = a' log m + &' , (8-25) 

As boundary conditions for the top of the chro- 
mospheres, we have chosen a velocity that is 
somewhat smaller than the sound velocity and 
a zero velocity for the minimum temperature. 
(The published models of Johnson do not con- 
tain microturbulence velocities.) 

For each photospheric model, we have cal- 
culated two linear chromospheres having dif- 
ferent gradients. One corresponds to a “hot” 
model (T o = 7000 K), and the other cor- 
responds to a “cold” model (T q = 4000 K). 
The boundary conditions are presented in Table 
8-1, where P is the total pressure. The exten- 
sion of the chromosphere, H chro , is also con- 
tained in the table. For illustration, Figure 8-3 
shows the behavior of the temperature versus 
column mass and the respective velocities. 

KE Line Formation of the Magnesium 
Resonance Lines. In order to study the main 
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Table 8-1 

Main Characteristics of the Chromospheres 


Effective 

Temperature 

Log g 

P q (dynes/cm 2 ) 

T 

o 

T . 

mm 

H (km) 

chro 

3500 

0.0 

icr 7 

7000,4000 

2380 

1.4 x 10 8 

3000 

2.0 

icr 5 

7000,4000 

2010 

9 x 10 5 



LOG M (G/cm 2 ) 


Figure 8-3. Distribution of temperature (solid lines) and velocities for the model 
3000/2.0/1.02. Velocities for the hot chromosphere (T q — 7000 K) are represented 
by crosses, and those for the cold chromosphere (T = 4000 K) by points. 


characteristics of the presented atmospheric 
models, we have calculated the KE line forma- 
tion of the resonance lines of Mg I at X2852 and 
Mg II at X2798. The latter line represents a 
mean of the h and k lines of Mg II. (The men- 
tioned IUE observations do not resolve the h 
and k lines.) The KE calculations have been 
made by means of the LINEAR code indicated 
earlier. This code has been adapted to the con- 
dition appropriate to carbon-rich atmospheres 
and has been modified to use Voigt profiles. 

One of the main characteristics of our KE 


calculations is that we assume the complete 
redistribution of photons. This analysis has 
been improved in studies considering partial 
redistribution. For our exploratory purposes, 
however, the assumption of the complete 
redistribution will be sufficient for the study of 
the core line formation. Only the indicated 
resonant lines are treated explicitly in the 
LINEAR code. All the other permitted transi- 
tions are considered as fixed. The radiation 
temperatures (T ) for the photoionizations 
and fixed transitions are equal to T ff (with a 
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dilution factor of 1/2) in the chromospheric 
layers where T is larger than T gff . In the 
following deeper layers, 7j. ad = 7\ A fixed- 
level atom is used for both Mg I and Mg II 
atoms. A solar Mg abundance is used. 

A strong absorption Mg I line at 2852.1 A 
is observed in many late-type stars. (This can 
be seen, for instance, in the spectrum of the M3 
III star, t) Gem, in Johnson and O’Brien 
(1983).) This line, however, appears to be ab- 
sent or weak in the spectra of the N carbon stars 
observed by Johnson and O’Brien. Normally, 
a strong absorption should be produced by the 
cooler photospheres of these stars, and this is 
not observed. Instead, a small emission feature 
appears at this wavelength in all the observed 
stars. Is this feature an indication of a “filling- 
in” process? Can an emission in the chromo- 
sphere produce this filling-in? All the Mg I lines 
obtained with these simple models are absorp- 
tion profiles which show no sign of a central 
emission feature because of the strong 
photoionizing control of the line-source 
function. 

An alternative filling-in process has been 
suggested by Johnson and O’Brien (1983) in 
which a large shell of gas could exist above this 
chromosphere. This shell must be sufficiently 
cool for Mg to be neutral and large enough to 
emit Mg I. The low-lying emission lines of Ti 
I, V I, and Zr I found by Gilra (1976) in the 
carbon star, UU Aur, could be formed in even 
larger shells because of their lower abundances. 
Scattering processes may also produce those 
emissions. 

Figure 8-4 shows the resulting theoretical 
profiles of the resonance line of Mg II at X2798 
for the 3000/2.0/1 .02 model. In both cases (hot 
and cold chromospheres), the line appears in 
emission with a central reversal. The shapes of 
these profiles are more or less typical of those 
found in late-type stars. One characteristic of 
this line is to be optically very thick, beginning 
to get thin only on the very top of the chro- 
mosphere. This line is, in all cases, formed at 
layers higher in the chromosphere than the 
Mg I line. The Mg II lines are much more con- 
trolled by collisions, indicating that very low 


radiation is available to photoionize from these 
two levels of the transition. The respective con- 
tinua are at X824 and XI 169, whereas those for 
Mg I are XI 621 and X3756. For a general discus- 
sion on collisional and photoionization control 
processes, see Jefferies (1968) and Athay 
(1972). 

The effects of the microturbulent velocities 
increasing outward are shown in Figure 8-5. 
The presented profiles are calculated for the 
same model 3000/2.0/1.02, both with and 
without microturbulence. The global effect of 
the velocity field is to broaden the line and shift 
the emission features, forming a broader rever- 
sal. More or less the same qualitative behavior 
has been found by Basri (1980) by means of a 
partial redistribution analysis. These turbulent 
effects are important in the detailed analyses 
of supergiants, which have much broader emis- 
sion cores than those of dwarfs or giants. 

The Avrett-Johnson Model. Apart from the Mg 

1 and Mg II characteristics discussed above, any 
chromospheric model calculation of a warm N 
carbon star should take into account other im- 
portant observational facts that act as con- 
straints in the modeling construction. These 
are: (1) the presence of the emission line of C 
II at 2325 A, (2) the notorious absence of the 
hydrogen Balmer lines (Yamashita, 1972, 
1975), and (3) the reversal in the Ca II H and 
K lines observed in R stars (Richer, 1975) but 
not in N stars. 

Avrett and Johnson (1984) recently con- 
structed semiempirical chromospheres to fit 
these observational constraints (with the excep- 
tion of the Ca II lines). The chromospheric 
model of an N star should be sufficiently cool 
to produce a low number density of the level 

2 of hydrogen (n 2 ) to prevent the production 
of an Ha feature and should, at the same time, 
be sufficiently hot to produce Mg II and C II 
line emissions. 

Their calculations are much more complete 
than the linear chromospheres calculated in the 
preceding subsections. They solved the 
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hydrostatic equilibria together with the KE of 
various elements such as H, H", H 2 , Mg I, Mg 
II, C I, and C II and considered partial 
redistribution in the calculation of the Mg II 
line. The radiation temperatures used in their 
calculations (between 2200 and 2500 K) are 
smaller than those of the linear models of the 
subsection Chromospheres of Carbon Stars 
which have lower electron densities. 

The resulting models 1 and 2 are presented 
in Figure 8-6, taken from the work of Avrett 
and Johnson (1984). The cool model 1 fills the 
condition of not producing the H feature; it 
seems, nevertheless, to be too cool to produce 
sufficient Mg II emission to fit the observed 
one. A negligible C II line is produced with this 
model. The hotter model 2 is also not conve- 
nient because it contains sufficient n 2 to pro- 


duce an emission or absorption Ha feature, 
depending on whether the source function is 
larger or smaller than B. 

In conclusion, both simple linear models 
and more sophisticated ones produce emission 
lines of Mg II. However, it is not clear whether 
future models will produce a C II emission 
feature or will maintain the H condition. 
(Avrett and Johnson (1984) are presently in- 
vestigating hotter lower chromospheric 
models.) Two possible states appear to solve 
this problem: a high carbon abundance or a 
hydrogen underabundance. 

Radiative Cooling and Energetic Balance. 

Several studies have recently appeared in the 
literature which treat the important problem of 
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and Johnson , 1984). 


the energetic balance in the atmospheres of the 
Sun and late-type stars (Linsky, 1980, 1981; 
Avrett, 1980). Two main processes exist in 
stellar atmospheres that maintain a temperature 
structure in equilibrium. These are the cooling 
and heating by radiative processes and the 
heating by nonradiative mechanisms. To obtain 
the net radiative contribution of each chemical 
element, the cooling and heating must be 
calculated for all the important bound-bound 
and bound-free transitions. In the 
chromosphere, a net radiative cooling flux 
results. This net cooling is a direct measure of 
the nonradiative (mechanical or magnetic) 
heating required to produce the chromosphere. 
In late-type giants, the lines of Ca II and Mg 
II are the main contribution to the net cooling, 
followed by H and H“. As we go from the Sun 
to cooler giants, there is a growing importance 
of Mg II as a net cooler agent in the higher 
layers of the chromosphere. 
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For a transition between levels u (upper) and 
/ (lower), the radiative cooling rate is (Vernaz- 
za et al., 1981): 

- *- 1". + «„ I* j, *> - 

(8-26) 

"A, I*. 

This expression can also be written using the 
net g radiative bracket (Thomas, 1960). Note 
that the same symbol, $, is used for two dif- 
ferent quantities: 

*'ul = hv n u A ul c * ( 8 ' 27 ) 

where 

S / dv 

6 = 1 . (8-28) 

s 1 

Figure 8-7 presents the values of $' u{ for the 
resonant line of Mg II calculated for the models 
3500/0.0/1.05 and 3000/2.0/1.02 with hot 
chromospheres. The cooling rates for the higher 
gravity model are found to be larger than those 
of the log g = 0.0 model. In fact, denser at- 
mospheres are better radiators than diluted 
ones. If we can express 4>' w/ as a function of 
the geometrical depth, z, the net cooling rate 
can be known by calculating the integral, 
j<F' (z) dz , for the positive values of $>' w/ . The 
cooling rates shown here are presented for il- 
lustration purposes. Realistic values can be ob- 
tained by matching theoretical and observa- 
tional integral fluxes. 

Ultraviolet Spectra Variability. Systematic 
observations of ultraviolet spectral time varia- 
bility are rare in S- and C-type cool giants (Cas- 
satella et aL, 1980, for the S-type star, X Cyg; 
Querci and Querci, 1985, for the N-type car- 
bon star, TW Hor). Interpretation of the car- 
bon star, TW Hor, appears to be very difficult 
because of the erratical behavior of emission 
lines such as those of Fe II at 3280 A. As is the 
case in these carbon stars, several constraints 
must be taken into account in the atmospheric 


model of TW Hor, such as the absence of Bal- 
mer lines (either absorption or emission) and 
the absence of emission peaks of Ca II H and 
K lines. The Mg II emission lines appear to have 
a more or less normal behavior, indicating 
smaller variability. (For a discussion on the 
observations of TW Hor, see M. Querci, this 
volume.) 

To understand the spectral characteristics of 
TW Hor, Querci and Querci (1985) analyzed 
a simple energetic budget which considers dif- 
ferent types of physical mechanisms such as 
short-period acoustic waves, variable magnetic 
fields, and flare-like activity. The emergent pic- 
ture that appears is that of a chromosphere with 
a flat temperature minimum and a gradual in- 
crease in higher layers. However, this increase 
of temperature is low enough to prevent emis- 
sion peaks in Ca II H and K lines, but strong 
enough in higher layers to produce ionic emis- 
sion features of the more abundant elements, 
Mg and Fe. Nevertheless, it is not clear how this 
chromospheric model will be able to prevent 
Balmer lines. The erratic behavior of the Fe II 
lines is, however, the most difficult aspect to 
understand. A variable La flux which induces 
variable Fe II emission lines (Furenlid, 1984) 
will not be an appropriate mechanism here be- 
cause of the atypical characteristics of H in 
these stars. Querci and Querci (1985) state that 
the cause of this variability can be due to varia- 
ble wave pulses (or variable magnetic field) in 
the chromosphere which is heated by short- 
period acoustic waves (Schmitz and Ulmsch- 
neider, 1981). The Mira S-type variable, X Cyg, 
observed by Cassatella et al. (1980) has been 
qualitatively interpreted using the dynamical 
shock-wave model of Hill and Willson (1979). 
(See the section Radial Pulsation Modes.) 

SHOCK- WAVE GAS DYNAMICS 
AND PULSATIONAL THEORY 

The preceding sections of this chapter 
discussed stellar photospheres and chromo- 
spheres as stationary components of static 
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Figure 8-7. Cooling radiative functions for the Mg II line . Only the positive values con- 
tribute to the cooling . The solid line represents the hot chromosphere of the log g = 0.0 
model. ( The corresponding scales are on the left and on the bottom of the figure.) The 
points represent the hot chromosphere of the log g = 2.0 model. (The corresponding scales 
are on the top and right of the figure.) The top of the chromosphere is on the left. 


atmospheres of single stars. However, some M 
and C late-type giants, as the long-period 
variables (Mira variables), present a complex 
time variability of their line spectra which re- 
quires an appropriate dynamical scenario for 
interpretation. Progressive shock waves are 
probably the principal physical mechanism that 
can explain this complex behavior. These 
shock-wave models are described by M. Quer- 
ci (this volume), who focuses mainly on the 
observational aspects. Here we will put more 
emphasis on some of the basic theoretical as- 
pects related to dynamical effects on the atmo- 
spheric structure and pulsational character- 
istics. 

First, we will present a short discussion on 
the principal physical basis of the classical 


theory of shock waves in order to analyze the 
current research on this field. 

Elementary Shock-Wave Theory 

Gas dynamics is characterized by quantities 
such as velocity (u), pressure ( P ), density (p), 
etc. which generally vary continuously. How- 
ever, discontinuities of these quantities, which 
will represent shock waves, can exist. These 
discontinuities can move with different veloci- 
ties of the gas flow itself and can be crossed by 
this flow of gas particles in their motion. If we 
consider a fixed coordinate system in which we 
define v n as the gas velocity normal to the sur- 
face and v s as the velocity of the discontinuity 
(shock wave), u = v n ~ v s will be the velocity 


394 




of the gas relative to the discontinuity. Let us 
call gas 1 the one into which the discontinuity 
moves (preshock) and 2 the gas which follows 
the discontinuity (postshock). The respective 
velocities will be w and u y 

One of the most important theoretical tools 
that enable us to treat shock waves is the fact 
that the laws of mass, momentum, and energy 
conservation can be applied to those discon- 
tinuity regions. Applying these laws (see, for 
instance, Zeldovich and Raizer, 1966), we have 
the following respective equations: 

G,«, = e 2 « 2 > (8-29) 


Since J 2 is positive, two possibilities may 
arise: P 2 >P V V, > V 2 or P 2 < P ] , K, <K r In 
realistic shock waves, only the first case may 
actually occur. This corresponds to stable com- 
pression waves that result in an increase of en- 
tropy (Zeldovich and Raizer, 1966). 

We can write Equation (8-31) in the follow- 
ing form: 


e + P. V + — J 2 V 2 
1 1 1 1 


= e. + P„ V. + — J 2 V 2 . 
2 2 2 -1 2 


(8-34) 


P, + Q^u] = P 2 + q 2 u\ , (8-30) 

P, U 2 p U 2 

e + — — + — — = e + — + — (8-31) 

e i 2 2 e 2 2 

where P, q , and e are the pressure, density, and 
internal energy, respectively. Equation (8-29) 
reflects the existence of a continuous mass flux 
crossing a “zero-mass” discontinuity. Equation 
(8-30) represents the existence of the continuous 
momentum flux, meaning that the force exerted 
by the gases on each other across the discon- 
tinuity must be equal. If no external energy 
sources are considered, Equation (8-31) means 
that the internal energy, e, of a given element 
is the result of the compressive work done on 
the element by the surrounding medium. 


Replacing J 2 from Equation (8-33), we find 
finally: 

e 2 - = \ (P,+P 2 )(V*V • ^ 

This relation is known as the shock adiabatic 
or Hugoniot relation. Graphically, this relation 
is represented in a PK-plane by Figure 8-8. An 
important conclusion can be obtained simply 
from this relation. First, we will determine the 
physical behavior of weak shock waves for 
small values of P 1 - P and V - V . Equation 
(8-33) can be written in first approximation as: 



The velocities u j and u being equal in the 
same approximation, we have: 


The Hugoniot Relations. A new important rela- 
tion can be obtained from the conservation 
laws. Introducing the specific volumes, V = 
1/q 1 and V 1 = 1 /g , we obtain from Equation 
(8-29): 

Q : u i = Q 2 u 2 = J ( 8 “ 32 ) 

that u { = JV X and u 2 = JV r Substituting in 
Equation (8-30), we obtain 



Wj = u 2 = u = JV . (8-37) 

With Equation (8-36), we obtain: 



which is the velocity of sound C for that 
medium. We can conclude that, in the first ap- 
proximation, the propagation velocity of a 
weak shock is equal to the velocity of sound, 
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Figure 8-8. Graphical representation of the 
Hugoniot curve in a pressure-specific volume 
plane. 

or in other words, the weak shock is a small 
perturbation that is approximately the same as 
a compression sound wave. Because real shock 
waves are compression waves, we must have 
P 2 P r The point, 1 (P { , Kj), must lie below 
the point, 2 (P 2 , V 2 ), in Figure 8-8. Chord 12 
has a slope J 2 larger than the tangent to the 
adiabatic at point 1: 


tangent at point 2, we obtain in the same man- 
ner as before that: 

u 2 <C 2 . (8-42) 

This means that the gas flows out of the 
discontinuity with a subsonic velocity, or in 
other words, the shock wave propagates at a 
subsonic velocity with respect to the compressed 
gas behind it. Because of the compression, q 2 
>e r we have by means of Equation (8-29) 
that Wj > u 2 . 

Shock Waves in a Perfect Gas. If we now con- 
sider the shock waves in a perfect gas with a 
constant ratio of specific heat, 7, we can in- 
troduce the internal energy: 

PV 

€ “ 7 , (8-43) 

7- 1 

into the Hugoniot relation (8-35) and obtain the 
following relation: 

— = 1 . (8-44) 

V i (7- 1 )-Pj +(t+ 1)^2 


> - 


dp 

dV 


(8-39) 


Multiplying both sides byK*, we have: 


dP 


(8-40) 


= C] 


Then, 


u, >C, . 


The gas flows into the discontinuity with a 
supersonic velocity. In other words, the shock 
wave propagates at a supersonic velocity with 
respect to the undisturbed gas 1 . The slope of 
the chord being smaller than the adiabatic 


Also, using the state equation of a perfect 
gas, we obtain the ratio of temperatures of the 
two sides of the discontinuity to be: 


T 2 f 2 ( 7+ l)/^ +P 2 (7-1)P 2 

T l ~ /> 1 (7-l)f J 1 +P 1 (7+l)i > 2 


• (8-45) 


From Equations (8-29), (8-30), and (8-44), 
we can obtain the velocities of propagation of 
the shock waves relative to the preshock and 
postshock gases: 


(8-41) 7 ^{(7- +(7+0^1- , 


(8-46) 


u\ = J ^-{(7+ l)/ > 1 +(7- ^ 
y /{(7-l)/ , 1 +(7+l)/ > 2 | 


(8-47) 
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In the case of very strong shock waves for 
which P 2 » Pj and P 2 » (7 + 1 )P { / 
(y - 1), we obtain the following result from 
Equations (8-44), (8-45), (8-46), and (8-47): 


"l 7+1 

K <?, 7-1 


(8-48) 


t _2 _ 7-1 ^ 
r, * 7+1 


", f'.l 


(8-49) 


(7+l)P 2 V, | V2 , (8-50) 


M 2 =[7 (7-l) 2 P 2 F 1 /(7+l)J 


1/2 


. (8-51) 


These relations have important physical im- 
plications. First, the ratio, 7^/7^, can be as 
large as P 2 /P x \ second, the velocities, u l and 
u v depend on the square root of P r This 
means that T 2 and the velocities can take very 
large values. 

This is not the case for the density, g r In 
fact, the ratio of densities S will attain a con- 
stant value, depending on the value of 7: S = 
4 for a monoatomic gas with 7 = 5/3, and S 
= 6 for a diatomic gas with 7 = 7/5. If full 
vibration/excitation is taken into account with 
7 = 9/7, 5 will be equal to 8. This ratio can 
be even larger if dissociation of molecules, ex- 
citations, and ionizations of the atoms or ions 
are taken into account. 

Shock Model Applications 


Three main directions characterize the work 
related to shock waves in the atmospheres of 
late-type stars: (1) the energy balance (the in- 
fluence on the thermal structure), (2) the modes 
of pulsation which induce the shocks, and (3) 
shocks as a mechanism that produces mass loss. 
As mentioned previously, the observational evi- 


dence of the presence of shock waves, their be- 
havior, and the models interpreting them are 
presented in detail in Chapter 2 (M. Querci). 
Here we will mention those observational as- 
pects related to specific points of the theoretical 
discussion. 

The motivation of the first theoretical ap- 
proaches was to explain the origin of the bright 
hydrogen lines. Gorbatskii (1961) made one of 
the first theoretical studies of the ionization and 
radiation of zones heated by shock waves. This 
study, however, did not take into account the 
dissociation of the H 2 molecule, which is an 
important cooling effect. Later, Whitney and 
Skalafuris (1963) studied the restricted case of 
the behavior of the high-temperature region be- 
tween 10000 to 30000 K near the shock front 
in a pure hydrogen atmosphere with no den- 
sity gradient in front of the shock. The effects 
of the dissociation of the H 2 molecule have 
been taken into account by Slutz (1976) in his 
study of the acceleration of a single isolated 
adiabatic (with no radiative energy losses) shock 
wave in a photosphere having a density grad- 
ient. This study shows that, if the shock has a 
large amplitude, the acceleration ends in a ter- 
minal velocity smaller than 25 km s" 1 (this 
value depending on the quantity of H 2 ). This 
could imply that, for stars with an escape 
velocity smaller than 20 km s -1 , an ejection of 
a shell of matter is possible. However, this kind 
of shock, called first shock, has been shown by 
Wood (1979) to be an atypical shock. In fact, 
a typical shock would be the one moving into 
matter falling inward, originated by a preceding 
shock, or matter flowing outward in a form of 
cool stellar wind. In addition, the proper initial 
conditions must be correctly established to con- 
sider the periodical nature of the problem. A 
sinusoidal varying pressure which simulates the 
star pulsation could be this initial condition. 

The intrinsic pulsation nature of cool giants 
or supergiant variable stars (see, for example, 
M. Querci, this volume) could be the source of 
shock waves. Even the microturbulent or mac- 
roturbulent velocities, in spite of being gen- 
erally ill-determined quantities, are sometimes 
found to be larger than the sonic velocities and 
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are therefore able to produce shock waves. 
However, it is in the interpretation of the com- 
plex spectra of Mira variables that the most 
powerful evidence of the presence of shock 
waves is found. The work of Willson (1976) is, 
in this respect, the most illuminating. Without 
considering in detail the complex behavior, 
doubling, and variations of absorption and 
emission lines (M. Querci, this volume), we 
mention the existence of two types of lines. The 
primary lines of Mg II, Ca II, Fe, and Ti II, 
formed in the region behind the shock (region 
2) and the secondary lines of Sc I, Mn I, and 
Fe I (fluorescent or pumped lines) formed in 
front of the shock (region 1). 

The main argument of Willson (1976) is that 
the fluorescent lines are formed by excitation 
of radiation emitted by the primary lines in a 
region of velocity discontinuity (shock wave). 
Each fluorescent line has its proper primary line 
exciting radiation approaching or receding with 
a characteristic velocity difference, producing 
an appropriate shift. Willson was able to show 
that, in general, the complex spectrum of 
typical Mira LPV atmosphere could be inter- 
preted by a velocity discontinuity of the order 
of 40 to 80 km s' 1 , thus providing strong 
evidence of the existence of a shock wave. 
Broadly speaking, the picture is the following: 
a moderately strong shock with a mean veloc- 
ity of near 50 km s' 1 appears at the base of the 
photosphere shortly after the visual light max- 
imum. As this shock travels outward, it ther- 
mally excites the primary lines behind the 
shock; these lines at their turn radiatively ex- 
cite (or pump) the fluorescent lines in the re- 
gions in front of the shock. 

It is important to realize that the fluorescent 
lines permit one to determine only a velocity 
discontinuity (Ac), not the actual shock-wave 
velocity ( v s ) with respect to the center of mass 
of the star. Positive or negative Ac indicates 
whether an exciting atom is approaching or re- 
ceding from the region in which the excited 
(pumped) lines are formed. In this way, Au 
0 will mean that the compressed region 2 is 
moving outward. In this region, as discussed 


in the section Elementary Shock- Wave Theory, 
u 2 = v 2 - v s , where v 2 is the gas velocity nor- 
mal to the discontinuity. From relation (8-42), 
we see that this difference is smaller than the 
speed of sound, C 2 , which is less than 10 km 
s" 1 for a gas at T< 10 4 K and a typical density 
of 10" 12 gr cm" 3 . Taking into account the 
limited precision of the measurements, we can 
say that v 2 = u s . In that condition, the veloci- 
ty u ] (velocity of the shock with respect to 
region 1) is a good measure of Au obtained by 
the fluorescent lines: «, = Au. If region 1 is 
falling toward the star, the actual shock velocity 
would be smaller than Au, v s = u 2 = Au - u,. 
On the contrary, if the material is flowing out- 
ward, the velocity of the compressed region will 
be larger, u 5 = u 2 = Au + u,. 

Influences on the Thermal Structure. The ex- 
istence of a supplementary kinetical energy due 
to shocks produces an extension of the entire 
scale of the atmosphere of a late-type giant or 
supergiant star. This extension can be a factor 
10 greater than the nonshocked atmosphere. 
For example, for a star with a radius equal to 
10 10 cm, this scale is similar to the stellar 
radius (Willson, 1976). The shock waves must 
then be considered as spherical shock waves 
which produce observable effects at greater 
distances from the stellar surface. In fact, 
shocks can persist for more than a period, and 
the product of the period and the shock velocity 
is of the order of the stellar radius. 

We also expect to find structural changes in 
the distribution of temperature and density. 
However, these changes are not clearly deter- 
mined at the present time because of the 
absence of a detailed calculation of the coupled 
interaction of the structure parameters and the 
detailed study of the shocks. A comparison be- 
tween plane-parallel and the more realistic 
spherical isothermal atmospheres can be found 
in Willson (1982). The main theoretical ap- 
proaches have been made using two different 
models: adiabatic (Wood, 1979) and isothermal 
(Willson and Pierce, 1982). The isothermal con- 
dition is expected to be satisfied in the lower 
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dense regions of the photosphere in which some 
absorption lines can be formed. Due to the 
higher density, any shock-heated material can 
cool more efficiently. This cooling is then pro- 
duced in a small distance compared to the en- 
tire scale of the atmosphere. The isothermal 
condition is less valid in the outer less dense 
regions. We then expect a progressive change 
toward the adiabatic condition. 

Using a pure adiabatic model, Wood (1979) 
found a resultant decreasing temperature dis- 
tribution corresponding to a mean heating of 
the order of 5000 K in relation to an initial 
hydrostatic atmosphere. However, this adia- 
batic model produces an unrealistic high mass- 
loss rate. Using the isothermal approximation, 
Willson and Pierce (1982) found, in an at- 
mosphere which suffers repeated shock pas- 
sages, that the minimum temperature attained 
between shocks is an increasing outward func- 
tion, contrary to the adiabatic case mentioned 
before. Due to the failure of the isothermal ap- 
proximation at densities lower than log q - -15 
(corresponding to a radius larger than 900 R © ), 
these authors can no longer describe the gas in 
the outer regions by a single temperature. In 
this picture we do not have a clear idea of the 
height at which the shocked model will attain 
a temperature low enough to allow the forma- 
tion of dust grains. When dust grains are 
formed, they will produce heating by gas/grain 
collisions. At the same time, opposite cooling 
processes will appear such as rotational cool- 
ing of molecules or a simple adiabatic expan- 
sion. A recent work by Tielens (1983) has con- 
sidered all these factors in a study of the flows 
of Mira variables driven by radiation pressure 
on dust grains. Tielens especially considered the 
direct coupling of the cooling of the gas with 
the velocity gradients in a two-component fluid 
composed by dust and gas interacting by colli- 
sions. This study is related to the existence of 
a stationary layer located between 5 and 10 R* 
from the star, in which matter can return to the 
star or be pushed outward. The existence of this 
layer at 800 K at 10 R, from which the maser 
features of SiO could originate, is proposed by 
Hinkle et al. (1982) and Hinkle (1983). 


Detailed Structure of Shock Waves. Some 
studies of the detailed physical parts of stellar 
shocks appeared recently in the literature (Gillet 
and Lafon, 1983, 1984, in a strictly theoretical 
point of view; Fox et al., 1984, combining 
theory and observations). (See also the subsec- 
tion Shock-Wave Gas Dynamics of Balmer 
Emission Lines.) A complete detailed analysis 
of the shock regions would be quite difficult, 
especially when the real finite shock front is 
considered, in which very large gradients of the 
physical variables have to be taken into ac- 
count. Appropriate mathematical methods will 
then be necessary. Perhaps one promising 
method for future studies is the “Adaptive- 
Mesh Radiation Hydrodynamics” (Winkler et 
al., 1984). Here, an adaptive coordinate system, 
fixed neither to the fluid nor to the laboratory, 
is free to follow the evolution of the flow and 
the radiation field. The importance of this 
method is that, although the observed quan- 
tities are measured, for instance, in a comov- 
ing frame of the fluid, the radiation transport 
equation can be written in this adaptive frame 
which can move with arbitrarily high velocities. 
Large gradients such as shock fronts can then 
be treated. 

As mentioned in the subsection Shock 
Waves in a Perfect Gas , the ratio S = q 2 /q { 
can change if excitation and ionization are 
taken into account. It can be demonstrated that 
the increase of S for a given T depends on 
molecular dissociation and the internal energy 
of the particles (e.g., ionization) (Zeldovich and 
Raizer, 1966). Values larger than S = 10 can 
be obtained for high degrees of dissociation and 
ionization. In a medium composed of a mix- 
ture of gases, we expect to find a variable 
behavior of 5 as the strength of the shock 
increases. 

In a medium formed by pure H, Gillet and 
Lafon (1984) analyze the precursor region (e.g., 
the immediate preshock region), putting em- 
phasis on the importance of molecules in the 
weakening of the shock. For this, they con- 
sidered the physical nonlinear relation between 
the precursor and the postshock region. In fact, 
before being affected by the advancing front, 
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the precursor is already influenced by the radia- 
tion coming from the hot postshock region. 
This radiation photodissociates the molecules 
and photoionizes the medium, producing a 
radiative cooling of the postshock region. This 
cooling is further accentuated by the thermal 
dissociation of molecules that cross the front. 

Pulsation Theory 

Radial pulsations are clearly seen in a large 
part of luminous M and C cool stars (see M. 
Querci, this volume). The aim of the theoretical 
research on the pulsation theory on these stars 
has been to answer the following fundamental 
and difficult questions: (1) Which is the pulsa- 
tion mode? (2) Is pulsation an effective mass- 
loss mechanism? To answer these questions, 
linear (and nonlinear) adiabatic (and non- 
adiabatic) pulsation models have been 
developed in recent years. However, even with 
the great effort of some authors, no definite 
answers to these questions have been given. 

Pulsations are produced in the large diluted 
envelopes of the stars in which the high 
luminosity determined by the core mass is be- 
ing transported essentially by a convective pro- 
cess. The luminosity in these late-type stars can 
be so high that a net positive energy produces 
a dynamical instability. The pulsation provokes 
shock waves that locally heat and probably ex- 
tend the atmosphere. This pulsational instabili- 
ty can relax into a steady regular pulsation such 
as those observed in LPV stars. However, the 
pulsations can also be transformed into strong 
relaxing oscillations in which mass loss can oc- 
cur (Tuchman et al., 1978). Nevertheless, the 
success of this mechanism in interpreting real 
variable stars unfortunately depends on a bad- 
ly known time-dependent convection theory. 
Another important theoretical input is the 
necessity to introduce realistic surface boundary 
conditions. Wood (1980) clearly shows these 
necessities when he discusses the uncertain 
results concerning the pulsational instabilities 
mentioned above and those of Wood (1979) 
and Willson and Hill (1979). In these models, 
not only convection is the weakest point, but 


also the mass loss produced. In fact, only few 
external zones are ejected, and the considera- 
tion of a proper external boundary, such as cir- 
cumstellar material, could inhibit this mass loss. 
A general discussion on these aspects can be 
found in Wood (1982). More details concern- 
ing mass loss can be seen in Goldberg (this 
volume). 

Nonradial Pulsations. Although pulsation has 
been discussed as a radial oscillation, there is 
a possibility of the existence of nonradial 
oscillations as supplementary energy- 
transporting mechanisms other than convec- 
tion. These oscillations have been explored by 
linear nonadiabatic calculations in late-type 
models by Ando (1976). Later, Smith (1980) 
considered the possibility that these nonradial 
oscillations could be the source of what is com- 
monly called “macroturbulence.” The pre- 
dicted amplitudes of this “late-type star 
seismology” are nevertheless very small, prob- 
ably a few tenths of ms -1 . They would 
therefore be very difficult to detect. Recent 
observations in this direction are not yet con- 
clusive (Smith, 1983, in a Tau and a Boo; 
Moon et al., 1983, in an M-S giant, e Oct). 
Other promising candidates could be the 
following M semiregular giants which show 
short-period variations, particularly of the spec- 
tral region around the Ca I X4227 absorption 
feature: R Crt (Livi and Bergmann, 1982) and 
L 2 Pup (Gomez Balboa and de la Reza, 1985). 

Radial Pulsation Modes. One of the fundamen- 
tal relations of radial pulsational stellar theory 
is the period-density relation (see for instance, 
Cox, 1980) which introduces the pulsation con- 
stant, Q , in the form: 

Q = P(q/qq) U2 

( 8 - 52 ) 

= P(M/Mq) w2 {R/RqY v1 . 

However, as has been shown in linear 
nonadiabatic models (Fox and Wood, 1982), 
Q is not a constant for any of the three lower 
order radial pulsation modes (fundamental, 
first and second overtone modes). Figure 8-9 
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Figure 8-9. Q values for the three lowest order 
radial pulsation modes plotted against period 
(from Wood, 1982). 

presents their corresponding values of Q q , Q v 
and Q 2 against the first overtone period. If 
molecular opacity is introduced, Q 0 will be 
reduced. On the contrary, Q } will be 
somewhat larger. Another result of Fox and 
Wood is that, for high-luminosity stars, the 
nonadiabatic fundamental period, P q9 is 
shorter than the adiabatic, P Q . This has an im- 
portant consequence because the adiabatic 
dynamical instabilities produced by large 
periods (Tuchman et al., 1978) disappear when 
nonadiabatic effects are included. In this linear 
nonadiabatic regime, the fundamental period, 
P Q , behaves as P Q a R 2 rather than P a R 3/2 of 
the classical relation (8-52). An approximate 
simple relation can be found for Q 0 between 
definite values of Mand P. No simple relation, 
however, can be found for the first overtone 
period, P , which behaves as R 3/2 for higher 
mass and as R 5 2 for lower mass. These results 
are found for typical galactic-disk long-period 
and semiregular variables. The L-M-R rela- 
tions change when massive disk supergiants (up 
to 25 Mq ) or other populations with different 
metallicities are considered. 

Linear and nonlinear nonadiabatic pulsation 
models have also been investigated by Ostlie et 
al. (1982) in view of applications to typical 
galactic Mira stars. They obtain P-M-R rela- 
tions for fundamental and first overtone pulsa- 


tions, examining in particular the effect of a 
turbulence pressure. In the linear case, pressure 
effects appear to be important for the fun- 
damental mode and lower mass. In the 
nonlinear models, turbulence amplifies the 
nonadiabatic effects. 

How do Mira stars pulsate? This is a con- 
troversial subject, and no definite answers have 
been given. Looking back to the past, Kamijo 
(1963a, 1963b) concluded that Miras do not 
pulsate in the fundamental mode, but rather in 
higher modes with shorter periods. Later, the 
research has been based on the possibility of 
finding stable pulsational adiabatic static 
models (Keeley, 1970; Wood, 1974, 1975). The 
introduction of more realistic nonadiabatic 
models has introduced stable models which can 
be applied to massive LPV supergiants and 
pulsating carbon giant stars. These stars show 
secondary periods with P ~ 2000 to 7000 days 
(Leung and Stothers, 1977), which are approx- 
imately 7 times the typical primary periods of 
300 to 700 days. Fox and Wood (1982) found 
several supergiant models in which P Q /P { ~ 7, 
indicating that the primary period is produced 
by excitation of the first overtone and the 
secondary period of the fundamental mode. 

Let us now discuss the arguments support- 
ing the fundamental mode pulsation of Miras. 
These are based on a dynamical model of shock 
waves (Willson and Hill, 1979; Hill and 
Willson, 1979; Willson et al., 1982). A general 
discussion can be found in Willson (1982). In 
particular, Hill and Willson (1979) developed 
a dynamical model, and in their calculations, 
they used the spherical symmetry and isother- 
mality assumptions. Their model is character- 
ized by the existence of two shocks. One lower 
shock emerges, with large amplitude, at the 
photosphere and lasts about a period. With the 
next period, another shock appears; two shocks 
can then be present: one stronger and lower and 
another weaker and higher. This lower shock 
is formed near the piston by the impact of the 
infalling large volume of rarefied preshock 
material with the rising material near the 
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photosphere. Because rarefied matter is falling 
and dense compressed material is rising, no net 
inward or outward motion is produced in this 
model. 

In this model, the lower shock gives rise to 
the hydrogen emission lines in the postshock 
region, whereas the passage of the weaker up- 
per shock model affects the absorption lines. 
The emission and absorption velocities then 
show the velocity variation across the lower and 
upper shock velocities, respectively. 

In both these shocks, we can define Av = 
v 2 - where v { and v 2 are the velocities de- 
fined in the section Elementary Shock-Wave 
Theory. Considering the relation (8-29), we 
obtain: 


-v 


2 


Av 

J^\ 


(8-53) 



(8-54) 


Because u 2 = v 2 - v $i where is the shock 
velocity, it can be written as: 


dr s 

V s = — = (1 +q)v 2 , (8-55) 


where 


Av/v 2 

q = . (8-56) 

S- 1 

For a strong shock (higher S), v s * v r 

One important consequence of the model of 
Hill and Willson (1979) is that the ratio defined 
by/3 = v 2 /v e is nearly constant, and equal to 
0.2. v e is the escape velocity (2 GM/r o ) w 2 , 
where r is the radial distance of the lower 

o 

shock. The value and behavior of /3 characterize 
the dynamical aspect of the model and is the 
basis of the interpretation of the observational 
results. In that case, 

v 2 (r) = pv e (r) = v 2 (r 0 ) (r/rj 1/2 . (8-57) 


For constant values of 0 and <7, the solutions 
of Equations (8-55) and (8-57) give: 


u 2 (t = P) = v, (t — 0) 


3 ”2 (°) p 

1 + - 0 +< 7 ) — “ 
2 r(o) 


- 1/3 


(8-58) 


Hill and Willson (1979) place, in an M-R 
diagram (see their Figures 10 and 1 1), the values 
of constant velocity v 2 (o) and v 2 (P) for the 
stars R Leo and O Cet. It is possible to find 
a region in this diagram that is best satisfied 
by the observational data. Probable masses and 
radii are then obtained for R Leo (M =1.8- 
4.0 Af © , R = 1.7 - 2.5 x 10 13 cm) and for 
Mira (A/ * 0.8 - 1.5 M© , R = 1.3 - 2.0 x 
10 13 cm). 

Lines of constant Q values also allow one 
to establish the best possible pulsation modes 
between Q = 0.1 (fundamental mode) and Q 
= 0.04 (overtone mode). Both values seem a 
priori satisfactory; however, considering sup- 
plementary data on masses and radii, Hill and 
Willson (1979) preferred the fundamental mode 
for both stars. 

The work of Hill and Willson (1979) has 
been criticized by Wood (1980, 1982). As men- 
tioned previously, considerable controversy ex- 
ists about this subject. Wood (1975, 1980) 
favors the first overtone mode. Wood derived 
larger radii (smaller Q) than those of Hill and 
Willson. Nevertheless, the principal arguments 
of Wood concern the dynamical and 
kinematical interpretations of the measured 
shock velocities in LPV. On the dynamical 
aspect, Wood (1979) found that, if a radial 
temperature gradient is introduced, the ratio 0 
is no longer constant, being 0.24 for a lower 
shock and 0. 1 1 for an upper shock. The other 
point of discussion concerns the possibility of 
measuring the values of the postshock velocity 
of the upper shock. Wood (1980, 1982) showed 
that, if the upper shock velocity can be mea- 
sured, its low accuracy will severely compro- 
mise the results. In fact, for a dynamical shock- 
wave model, Q is very sensitive to the 
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velocity values because Q a (u |owei A upper ) 3 , and 
an uncertainty of ± 1 km s' 1 will introduce an 
error in the determination of Q of the order of 
± 0 . 1 . 

In conclusion, we can say that appropriate 
surface boundary conditions will help future 
studies. These conditions will probably affect 
primarily the overtone pulsations, which are 
more surface phenomena than those of the fun- 
damental mode. A new approach, based on the 
detailed shock evolution of the Balmer emis- 
sion lines, is now beginning to be explored. 


Shock-Wave Gas Dynamics of Balmer Emis- 
sion Lines. The observed behavior of the 
Balmer emission lines in different phases of 
Mira variables has been known since the 
pioneer works of Merrill (1945) and Joy (1947), 
but not until now have careful high spectral 
resolution observations of these lines been made 
in order to study the pulsation properties (Fox 
et al., 1984). A qualitative study of the dynami- 
cal evolution of Ha in the star, o Cet, has also 
appeared recently (Gillet et al., 1983), in which 
the variations of the Ha emission profile is in- 
terpreted as a spherical shock wave propagating 
outward from the photosphere. The shock 
leaves the photosphere when Ha emission ap- 
pears for the first time in the cycle. Because of 
the occultation by the photosphere, the 
observer sees only the advancing part (detec- 
ting a blue-shifted emission). The emission pro- 
file appears to be mutilated by absorptions of 
an external molecular region. This mutilation 
diminishes as the shock leaves the molecular 
region. At that moment, the shock is quite far 
from the photosphere, and the observer is able 
to see the receding part of the always-outgoing 
shock, now producing a red-shifted emission. 
The atmosphere is then crossed from the 
luminosity maximum to the minimum by a 
single strong spherical shock wave which pro- 
duces the emission lines in the hot postshock 
region. However, observed evidence of a se- 
cond supersonic infalling shock in the same star 
has been detected by Ferlet and Gillet (1984) 


by means of an inverse P Cygni profile of Ti 
I. This infalling shock could be the return of 
nonejected upward-accelerated matter. 

Fox et al. (1984) published a series of high- 
resolution spectral observations of the Balmer 
lines Hy, H<5, Hf, and Hr? of nine Mira stars 
during half cycle around maximum phase. 
Observations of the thermal excited line, Mg 
I X4571, and the pumped line, Fe I X4202 
(pumped by the Mg II resonance line at X2795), 
were also made. The Balmer emission lines pre- 
sent a full base width of near 80 km s' 1 which 
decreases at later phases. H5 seems to be the 
most indicated line for the study of evolution 
of emission regions because it is the least af- 
fected by superposed absorptions. 

A detailed shock model to be applied to 
these observations is in progress (Fox and 
Wood, 1985). However, in advance of the 
presentation of this detailed model, a schematic 
simple shock is presented by Fox et al., (1984). 
This shock is shown in Figure 8-10. The precur- 
sor is formed by H and He at T ~ 2500 K, with 
an infalling velocity of 20 km s' 1 . The hot 
postshock region is at T = 34500 K. It is in this 
region that the broad emission lines are formed 
by an important Balmer line scattering. In fact, 
the broadening temperature of the lines indicate 
temperatures of the order of 30000 K. The 
velocities shown in Figure 8-10 agree with those 
obtained from high-infrared spectroscopy of 
CO lines (Hinkle et al., 1984). In the shock 
model of Fox and Wood (1985), the far side 
of the outgoing shock does not contribute to 
the observed Balmer emission because of the 
total absorption of the inward Balmer radia- 
tion. This is contrary to the Ha interpretation 
of Gillet et al. (1983) mentioned previously. The 
metallic lines observed are very important in 
that they permit the Balmer emission variations 
in the cycle to be related to the velocities ob- 
tained from these metallic lines. This relation 
suggests a slowing down of the shock with 
phase. A further study is necessary to determine 
where the Mg I and Fe I lines mentioned are 
formed. Are they formed in the precursor or 
in the postshock region? The detailed model 
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Figure 8-10. A schematic diagram of shock structure in a typical Mira variable near- 
maximum light. Velocities are in the rest frame of the star (from Fox et al., 1984). 


of Fox and Wood will possibly give the answers 
and may give new information concerning the 
pulsation modes of the Mira stars investigated. 
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